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ABSTRACT 



Methods, systems, and apparatus for mapping a nonspectral 
representation of a target color, such as an input color tuple, 
to a set of concentration values for a set of device-specific 
colorants. The invention includes using the input color tuple 
to derive a first set of colorant concentration values from a 
color lookup table; and refining the first set of colorant 
concentration values by an iterative non-linear process to 
generate a final set of colorant concentration values. The first 
set of colorant concentration values can derived by using an 
input color tuple as an index to obtain grid-point concen- 
tration values at two grid points of the color lookup table and 
calculating the first set of colorant concentration values as a 
linear interpolation of the grid-point concentration values. In 
an implementation that provides a color function table of 
color-additive function values, the iterative non-linear pro- 
cess can include iteratively (a) calculating an interim color 
tuple from an interim set of colorant concentration values 
and the color function table, the initial interim set of colorant 
concentration values being the first set of colorant concen- 
tration values, and (b) deriving an interim set of colorant 
concentration values from a difference between the input 
color tuple and the interim color tuple. The calculations can 
include calculating a partial derivative of an error function 
from the difference between the input color tuple and the 
interim color tuple; and using the partial derivative to derive 
a successor interim set of colorant concentration values from 
a current interim set of colorant concentration values. 

31 Claims, 3 Drawing Sheets 
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SPECTRAL COLOR MATCHING TO A 
DEVICE-INDEPENDENT COLOR VALUE 



BACKGROUND OF THE INVENTION 

The present invention relates to the fields of color match- 
ing and color management systems. 

For the rendering of color images from a source in which 



color tuple and the interim color tuple. The calculations 
includes calculating a partial derivative of an error function 
from the difference between the input color tuple and the 
interim color tuple; and using the partial derivative to derive 
a successor interim set of colorant concentration values from 
a current interim set of colorant concentration values. 

In general, in another aspect, the invention provides 
methods, systems, and apparatus, including computer pro- 
gram apparatus, for color matching. The invention includes 



the color information is in a device-independent format, yet JQ a p para tus or steps that receive a target color value, the target 



to be rendered in a device-specific manner, conventional 
device -independent color management uses calorimetric 
information only. Device-independent triples, such as XYZ, 
are converted to colorant concentrations for printing on a 
particular device using a lookup table indexed by a device - 
independent triple whose lattice points contain colorant 
concentrations. Because lookup tables used for this purpose 
must necessarily be of a practical size, image pixel values 
generally miss the lattice points. To compensate for this, a 
linear interpolation is conventionally carried out to approxi- 
mate the intermediate colorant values. The lookup table 
lattice points must have very accurate values for both the 
primary colorants and their corresponding overprints. Table 
sizes in the ranges 10x10x10 to 33x33x33 are in common 



color value being a nonspectral representation of a target 
color; receive a spectral definition of each colorant in a set 
of colorants; obtain an initial set of colorant concentrations 
for the set of colorants, each concentration corresponding to 
15 a respective colorant in a set of colorants; and generate 
iteratively a final set of colorant concentrations by calculat- 
ing a rendered color value, the rendered color value being a 
nonspectral representation of a rendered color produced by 
a mixture of colorants defined by the set of colorant 
20 concentrations, calculating a difference between the ren- 
dered color and the target color, and modifying the set of 
colorant concentrations according to the calculated differ- 
ence. 



Advantageous implementations of the invention can 

use. In order to populate such tables, thousands of accurate 1$ fa^fa one 0 r more of the following features. A threshold is 

colorimetric measurements must be taken. received, and the iteration terminates when a calculated 

A number of difficulties are in the use of conventional difference between the rendered color and the target color is 

color management techniques. It is difficult to calculate less than the threshold or an iteration limit is reached, 

colorant concentrations that closely and perceptually match Obtaining an initial set of colorant concentrations includes 

a requested device -independent value. Because assumptions 30 using the target color value to select an index into a colorant 

about the linearity of the given color space are usually concentration lookup table; and obtaining the initial set of 

incorrect, it is difficult to obtain accurate interpolation. colorant concentrations from the indexed contents of the 

Extensive measurements must be taken to fill a lookup table. lookup table. Calculating a rendered color corresponding to 

It is difficult to produce invariant matches, that is, colors that the set of colorant concentrations includes calculating a 

will match under all lighting conditions. It may be necessary 35 rendered color reflectivity function from the set of colorant 



40 



concentrations and a set of colorant reflectivity functions, 
each colorant reflectivity function corresponding to a respec- 
tive one of the colorants; and calculating a rendered color 
value from the rendered color reflectivity function. The 
rendered color value is calculated using an illuminant 
selected in response to a user input identifying desired 
viewing conditions. Modifying the set of colorant concen- 
tration values includes calculating a set of colorant concen- 
tration corrections from the difference between the rendered 



to limit color management to CMYK device spaces and to 
require that many overprints be measured to populate the 
lookup tables. 

SUMMARY OF THE INVENTION 

In general, in one aspect, the invention provides methods, 
systems, and apparatus, including computer program 
apparatus, for mapping a nonspectral representation of an 
input color to a set of concentration values for a set of 

device -specific colorants, where the representation is an 45 color value and the target color value. The target color value 

input color tuple. A tuple is a representation of two or more is a triple defining a location in a color space. The target 

components, such as XYZ color components. The invention color value defines a device-independent color. The difler- 

includes using the input color tuple to derive a first set of ence between the rendered color and the target color is 

colorant concentration values from a color lookup table; and calculated in a perceptually linear color space. The percep- 

refining the first set of colorant concentration values by an 50 tually linear color space is CIELAB or CIELUV. The set of 

iterative non-linear process to generate a final set of colorant colorants is an arbitrary set of two or more colorants. The 

concentration values. spectral representation of each colorant derives from reflcc- 

Advantageous implementations of the invention include tivity measurements of the colorant taken at intervals over 

one or more of the following features. The first set of the visible spectrum without overprinting. Modifying the set 

colorant concentration values is derived by using an input 55 of colorant concentration values includes modeling the 



color tuple as an index to obtain grid -point concentration 
values at two grid points of the color lookup table and 
calculating the first set of colorant concentration values as a 
linear interpolation of the grid-point concentration values. In 
an implementation that provides a color function table of 60 
color-additive function values, the iterative non-linear pro- 
cess includes iteratively (a) calculating an interim color 
tuple from an interim set of colorant concentration values 



difference between the target color and the rendered color as 
a scalar field over the set of colorants; and calculating a 
partial derivative of the surface with respect to each colorant 
to calculate a correction for each colorant. 

Advantages that can be seen in implementations of the 
invention include one or more of the following. In many 
instances, use of the invention achieves an invariant, that is, 
a spectral, match to target color. Successful operation of the 
invention to converge to a matching color does not depend 



and the color function table, the initial interim set of colorant 

concentration values being the first set of colorant concen- 65 on interpolation or colorimetry. If a color lookup table is 

tration values, and (b) deriving an interim set of colorant used to accelerate convergence to a color match, it may be 

concentration values from a difference between the input sparse and approximate. The invention can produce colorant 
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concentrations that closely and perceptually match a space. In alternative implementations, target color value 110 

requested device-independent value. The invention works is represented in LAB or any of the other device- 

with any transparent colorants and any number of them. The independent color spaces derived from the XYZ color space, 

invention does not require a spectrum for the target color — it or the color value 110 is represented in an form that can be 

only requires calorimetric data for the match— and never- 5 mapped to Lab or XYZ color space, 

the less it can generate a spectrum for the target color. The Colorant lookup table 111 can be implemented as an 

invention relates color management and color matching to indexed list of sets of colorant concentrations. Each element 

each other in a new way. °f tne ^ l contains a set of colorant concentrations that 

rr» j * 'i c u j- • specifies a mixture with a color matching the index color. In 

Iht , details of ^one or more embodiments of the , invention 0 F ne imp i ementationj lhe coloranl table 111 is 

are set forth in the accompanying drawings and the descnp- 30 ^^J^ in the color space in which the target 

tion below. Other features and advantages of the invention color vahie no is ^presented. In an alternative 

will become apparent from the description, the drawings, implementation, the table is always indexed by L*a*b* 

and the claims. coordinates. The colorant lookup table 111 returns a set of 

colorant concentrations. 

Colorant characterization data 112 embodies a set of 

FIG. 1 is a block diagram of a system for obtaining a measurements of reflectance of visible light off of each of 

spectral color match for a target color using any number of me colorants in a colorant set. The colorant characterization 

colorants. data defines the colorant set to the color matching engine 

n^- a u * c ™ *u a* u* • ~ ♦ i i 100. The colorant set can be any arbitrary set of two or more 

FIG. 2 is a flowchart of a method to ob tarn a spectral color 9n _ . / ., J . . A . 

. u * * ♦ i uf i. colorants. It is advantageous to provide the colorant char- 
match for a target color using any number or colorants. . . , ^ * r , j j ■ • 

° actenzation data 112 in the form of a color-additive function, 

FIG. 3 is a flowchart of a method to obtain the reflectivity which ^ ^ described beloWj for each ^lonM in lhe 

for a mixture of colorants, given the reflectivity of each of colorant set 

the colorants, j n 0Qe implementation reflectances are measured in 10 

FIG. 4 is a schematic diagram illustrating color manage- 25 nm mcrements ranging from 400 nm to 700 nm wavelength 

ment workflow in accordance with the invention. and at 25%, 50%, 75% and 100% of colorant down. In other 

Like reference numbers and designations in the various implementations, reflectances are measured over a different 

drawings indicate like elements. range of wavelengths, or using a different wavelength 

increment, or with a nonuniform wavelength interval. Colo- 

DETAILED DESCRIPTION 30 fant re fl ect i v ity data can also be measured for different 

FIG. 1 is a block diagram depicting an implementation of percentages of colorant down, ranging from just above 0% 

a system for obtaining a spectral color match for a target colorant down to 100% colorant down. The increment can 

color using a set of any number of transparent colorants. For be differ from 25% intervals, and need not be uniform, 

example, the set of colorants can be a set of CMYK inks or 35 Viewing conditions 113 specifies the illumination under 

a set of Pantone® Hexachrome® inks. (Hcxachrome is a which the rendered color is to be viewed. The viewing 

six-color printing process developed by Pantone, Inc. of conditions 113 can be default viewing conditions or can be 

Carlstadt, N.J.) Other colorants can also be used; and the specified by user input. The viewing conditions can desig- 

method for matching the target color by a mixture of nate one of the CIE-defined illumination functions, such as 

colorants is not limited by, nor dependent upon, the number 4Q D65. In the present implementation, the viewing conditions 

of colorants or the colors of the colorants used. include as observer the 2-degree standard observer. Other 

As shown in FIG. 1, the system has a color matching viewing conditions can also be selected, 

engine 100 receiving as input a target color value 110, an The threshold 114 specifies a maximum acceptable AE, 

optional colorant lookup table 111, colorant characterization that is, a maximum acceptable difference between the target 

data 112, an optional viewing conditions 113, and a thresh- 45 color value 110 and the rendered color produced by a 

old (or target error value) 114. The color matching engine mixture of colorants mixed according to the calculated set of 

100 first calculates and then outputs a set of colorant colorant concentrations. The threshold 114 can be a default 

concentrations 120 that defines a mixture of colorants that threshold value or it can be specified by user input, 

produces a color that matches the target color specified by The set of colorant concentrations 120 is generated by the 

target color value 110 under viewing conditions 113. 50 color matching engine 100 to specify the composition of a 

Optionally, the color matching engine 100 also outputs a mixture of colorants with a color that matches the target 

quality of color match 121, which is a measure of the degree color represented by the target color value 110. 

to which the rendered color matches the target color. Here, The quality of color match 121 is generated by the color 

and throughout this specification, the term rendered color is matching engine 100 to indicate the closeness of the match 

used to refer to the color produced by a mixture of colorants 5S between the rendered color and the target color represented 

mixed according to a set of colorant concentrations. by the target color value 110. 

Optionally, the color matching engine 100 also outputs color The rendered color reflectivity 122 is a calculated set of 

reflectivity data 122 for the rendered color. reflectances of visible light off of the rendered color in 10 

Target color value 110 is a nonspectral representation of nm increments ranging from 400 nm to 700 nm wave- 

the target color to be produced by the color matching engine 60 lengths. In other implementations, different wavelength 

100. A nonspectral representation is one that does not ranges as well as different wavelength increments can be 

include spectral data such as reflectivity data. For example, returned, and the wavelength interval between reflectances 

an representation of a color as tri-stimulus values is a can be nonuniform. To avoid having to interpolate, mea- 

nonspectral representation. In the. implementation that will surement should be taken at the wavelengths for which the 

be described in detail below, target color value 110 is 55 viewing conditions power spectrum P x vector has values, 

represented in as values (X, Y, Z) of the 1931 CI E As shown in FIG. 2, the color matching engine imple- 

(Commission Internationale de l'Eclairage) XYZ color ments an iterative process to generate a set of colorant 
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concentrations for a mixture of colorants to match target 205 and 206), the color matching engine outputs the gen- 
color value. The color matching engine receives as input a erated set of colorant concentrations and, optionally, the 
target color value (step 200). The color matching engine quality of color match and the rendered color reflectivity 
determines an initial set of colorant concentrations (step (step 209). The quality of color match is represented by the 
201). The initial set of colorant concentrations can be 5 difference calculated at step 204. The rendered color reflec- 
determined as the indexed contents of a colorant lookup tivit y 122 is me reflectivity of the colorant mixture calcu- 
lable, where the index is selected based on the target color lated at ste P 202 for a mixture of colorants mixed according 
value. The initial set of colorant concentrations can be t0 the generated colorant concentrations, 
determined in other ways as well, such as by solving the FIG * 3 illustrates a process for calculating the reflectivity 
matrix equation (3), discussed below, in a least squares 10 of a colorant mixture 10 implement slep 202 (FIG. 2). Step 

202 requires as input colorant characterization data and a set 



sense. 



25 2JP A 



~, r *u i . i_- . iU . 4 4 . of colorant concentrations. At step 301, the color matching 

The core of the color matching process is the iteration rtU * • i * u ♦ • a * c u c.u 

, . . 4 , , .. *v 4 „ no „ . engine obtains colorant characterization data for each of the 

beginning at step 202 and ending at step 208. For step 202, f . . e - . .... - . - , - 

*u i * u- inA 5 . . 4 , c colorants in the form of a color-additive function for each of 

the color matching engine 100 requires as input a set of ^ co i orants 

colorant concentrations. Initially, the set of colorant concen- 15 « . • , . , . „ 

, . , . «• Color-additive functions are denved from colorant reflec- 

trations is determmed at step 201; m subsequent iterations, j * j- * *u i * * ^ u n i 

it _ 4 * , . ■ j j. * j tivity data according to the single-constant Kubelka-Munk 

the set of colorant concentrations is determined at step .208. . « ' , , ~ . ° ..... Jr n i « 

™. . . . , , « • • . r model. The color-additive functions 8 X (which are K/S 

The color matching engme calculates reflectivity data for a . v . „ * M . 4l _ x i i * j r *i_ * 

ri T . , J ti i- . *• m Kubelka-Munk theory) are calculated from the colorant 

mixture or colorants mixed according to the concentrations a . . , . g . , £ . 

.„ JLI . * 1 • / in reflectivity data (a set of spectral measurements of the 

specified by the given set of colorant concentrations (step 20 , . ^ j- * *u c n * v u n w i • i 

™ J . f. ... . i j * • i i.j- colorants) according to the following Kubelka-Munk ink 

202). The method by which reflectively data is calculated is eaU ation- 
described more fully below, in reference to FIG. 3. 

The color matching engine then calculates a rendered 
color value representing the color rendered by a mixture of 
colorants mixed according to a set of colorant concentrations 
(step 203). The rendered color value colorant vector is 

calculated in L*a*b* color space (CIELAB) and is repre- where 

sented by coordinates (L', A', B 1 ) in that space. The color denotes the i-th colorant (note that the superscript does 

matching engine then calculates the difference between the not denote an exponent), 

vector (L, A, B) representing target color value in L*a*b* 30 ^ denotes a wavelength, and R x ' denotes the reflectivity of 

color space and the vector (L', A B') representing the the iAh co l 0 rant at wavelength X. 

rendered color produced by the mixture of colorants (step ^ an impi^u parameter in e x \ namely the percent 

204). The color difference is generally calculated as: color down at wh i CD the corresponding R^' were measured. 

Ag°^ (L-L') 2 +(A-A') z +(B-B') z (1) 35 10 ° tDer WOrdS ' there ^ 0De VeCt0r ° f lheta values for eacn 

colorant down percentage; in the implementation being 

As will be seen, it is advantageous to work with the square described, these are 25%, 50%, 75%, and 100%. 

of the color difference, calculated as: An interesting property of these functions 6 X * is that they 

are color-additive and thus the result of adding given con- 

&e-(l-L') 2 +(s\-a*) 2 +(b-b¥ (2) centrations of the colorant functions yields a colorant tunc- 

Other measures or color spaces, such as CIELUV, can be ,ion for j he , r 1 esultant c ° lon , Fo ' «* m P'»' ^Junctions for 

used. It is advantageous that the color space and color ^ an . and f 11 ™ ma y * calculated and added to yield the 

distance metric be perceptually uniform, at least in the range ^f 00 f ° r ,he 1 resul,ant hu % °f S reen - In the present 

of applications in which the color matching engine will be notatlon - th,s m y be wntten as follows - 

used. 45 

The color matching engine then determines whether the c '^ 1 

calculated difference is less than the threshold (step 205). If iGiCiYl 
the calculated difference is not less than the threshold, the 

color matching engine checks whether the calculation where { ranges over the ^ co^^ C ( cyan ) ^ ^ 

should be terminated (step 206). The calculation can be 50 (yellow) and A ranges over the visible spectrum, and where 

terminated, for example, if the error AE does not improve c . ^ the concen tration of the i-th colorant being used in the 

after 3 iterations, or after some maximum number of itera- mixture , in this case, the concentration of cyan or yellow, 

tions has been completed. More geae rally, where there are n colorants in a set of 

If the calculation is not terminated ("no" branch of colorants P (a set of primaries) and the colorants are denoted 

decision step 206), the color matching engine calculates for 55 by numbers i through n, the color-additive function 6 x m of 

each of the colorants the gradient of the color difference with a mixture m is given 5y lhe fol i owing equat io n . 
respect to each colorant concentration (step 207). The color 

difference gradient is then used to calculate a new colorant ^ 

concentration for each colorant in the set of colorants (step ^ = 2j Ci&x 

208) and the calculation continues at step 202. The gradient 60 
and colorant concentration calculations (steps 207 and 208) 

will be described in detail later. In performing this calculation, which vector of theta 

The iteration ends when the color difference is less than values 0 X ' is chosen depends on the value c ; . In the present 

the threshold ("yes" branch of decision step 205) or the implementation, the vector whose percent colorant down is 

calculation is terminated for some other reason, such as a 65 closest to c ( - is chosen. In an alternative implementation, the 

failure to converge ("yes" branch of decision step 206). vector values are calculated by interpolation for each colo- 

When the process ends ("yes" branches of decision steps rant according to c ; the concentration of the colorant. No 



06/29/2004, EAST Version: 1.4.1 



US 6,646,763 Bl 



particular kind of interpolation is required. For example, a 
linear interpolation or an interpolation along the gamma 
curve of the colorant in question can be used. 

The color matching engine obtains a set of colorant 
concentrations (a vector of c ( ) at step 303. Initially, the 
vector is provided by step 202; in subsequent iterations, the 
vector is provided by step 208. The color matching engine 
calculates a color- additive function 0 x m for a mixture of the 
colorants (step 302) by multiplying the vector of c ( - by a 



8 

-continued 



where k is a constant calculated so that bright white has a Y 
value of 100. 

Finally, the color matching engine calculates the mix- 
ture's coordinates (L\ A', B') by applying the well-known 



matrix whose columns consist of the colorant color-additive 10 linear transform to convert form XYZ color space to L*a*b* 



functions obtained in step 301. Mathematically, the color 
matching engine calculates a color-additive function for a 
mixture m of colorants according to the equation 
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(3) 15 



RT = I + S- y/s 2 + 2S 



where 



color space. 

Steps 207 and 208 of FIG. 2 will now be described in 
more detail. 

The gradient of the color difference squared, AE 2 , with 
respect to the concentration of the i-th colorant is calculated 
as: 



where the columns of the matrix are color-additive functions 
of the component colors derived from spectral data taken at 
10 nm intervals as described above. (Recall that the columns 
1 through n for the n colorants will be chosen to correspond 
in percent colorant down to the c,.) This is generally an 
overdetermined matrix equation, and it can be solved by 
known numerical methods, such as the method of normal 
equations, to achieve a solution in the least squares sense. 
However, the results achieved by following this approach 
are not satisfactory, and in the present implementation, the 
vector of c f values is calculated by iterative refinement — in 
particular, by searching for a minimum of the error term 
AE 2 , as will be described — rather than by using traditional 
numerical methods. 

The color matching engine then inverts the color-additive 
function of the mixture to obtain reflectivity data for the 
color rendered by the mixture (step 304). This is done by 
solving the Kubelka-Munk color-additive function for the 
reflectivity, and is expressed as: 



dA£ 2 



(4) 



20 



where 



W + 0* 



25 



30 



35 



■20? 



1000(4 - A') 
400(5-^) 

232(L- V) - 1000(A -A') + 400(1 - g) 



45 



L, A, B are the L*a*b* coordinates of the target color 
L', A', B* are the L*a*b* coordinates of the mixture color 
X n , Y n , Z n are the XYZ coordinates of the viewing 

conditions white point and 
X, Y, Z are the XYZ coordinates of the mixture color. 
The color matching engine uses the gradient of the color 
difference with respect to a colorant concentration to calcu- 
late a new colorant concentration for each of the colorants 
(step 208). More particularly, for each colorant, the value 

\ d&E 2 \ 



A process will now be described that can be performed at 
step 203 (FIG. 2) to calculate the coordinates (L, A, B) of a 
color rendered by a mixture of colorants mixed according to 
a set of colorant concentrations. Step 203 requires as input 
a colorant mixture's reflectivity data (R^" 1 ), tne CIE percep- 
tual color-matching functions (x x , y^, z^), and an illuminant 
P x (the power distribution of the viewing conditions white 
point). 

The color matching engine calculates the mixture's coor- 
dinates in XYZ space using discrete analogs of the well- 
known CIE integrals, namely: 



50 is scaled if a large mismatch occurs — for example, if a 
correction to a concentration has a magnitude greater than 
1.0 or if the resulting concentration cannot be realized 
physically. The scaled value is calculated to be no greater 
than 0.5. The concentration c,- is increased or decreased 

55 

d&E 2 ) 

according to the sign of — — 



60 by the original value or, if one has been calculated, by the 
scaled value. Because use of a lookup table provides good 
initial colorant concentration estimates, it is sufficient for 
generating the scaled value simply to divide the magnitude 
by 10 until it is less than one (if greater than one) and then 

65 divide by two until less than or equal to 0.5 (if greater than 
0.5), and to treat each of the scaled values as independent of 
the others. 
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The expression in equation (4), above, for the partial 
derivative of AE 2 was derived from the following formulas 
that relate the XYZ to L*a*b coordinates. 

A-a*»500(CAr/Y„) 1/3 -(r/r„) 1/3 ) 

fl-6*-20o((y/y„) 1/3 -(z/z „) 

While satisfactory, use of these formulas may result in less 
than optimal matches in areas of highlights, or for very light 
colors generally. To generate optimal matches for such 
colors, the partial derivative of the error function used for 
light colors is derived from the following formulas. 

i-x-'-iiecy/yj^-ie for r/r„>o.oos856 

L-L*«903.3(y/7„) otherwise 
A-a*-5Q0(f(X/X„)-f(YfY n )) 
B-b*-2Q0WYfY n )-f{ZfZj) 



fttyt 1 ' 3 for />0.0O8856 
/(f)-7.787f+l 6/116 otherwise 

The partial derivative for very light color has the same form 
as shown in Equation (4), above, with the following different 
values for the following terms. 



AA = 



1000- /-(A -A) 



AS 



AL 



400-/(5-^) 



232- 1 ■ {L-L l )~ 10O0- h {A- A 1 ) +400 : l : iB-^) 



where 



1=7.787. 

As will be described in more detail in reference to FIG. 4, 
one advantage of the present invention is that it can be 
implemented to provide a workflow that is backward com- 
patible in operation with the conventional workflows of 
widely used publishing systems. In a conventional 
workflow, a device independent triple is converted to the 



vided in that interpolated values can be provided as an 
output of the process as well as being used as a starting point 
for the non-linear process of iterative refinement 450 that 
has been described. This refinement process 450 can be 
implemented using a second precalculated table 440 — the 
table of color-additive function values — to refine the values 
produced in the first stage from the color lookup table, as has 
already been described. Briefly, the input color 420 and 
viewing conditions 405 are processed (operation 410) to 
generate an L*a*b* representation of the requested or target 
color 415. The color-additive function values 445 from the 
precalculated table 440 and a colorant vector C (either the 
initial vector 435 or a refined vector 480, which will be 
described) are processed to compute reflectances 460 of the 
15 colorant mix represented by the vector C (operation 455). 
The color represented by reflectances 460 and the requested 
color 415 are compared (operation 465) and the difference 
components 470 are used to calculate a new colorant vector 
C 480 using the techniques described in reference to Equa- 
20 tion 4, above (operation 475). If the vector C is close enough 
to the requested color 415, it is an output result 485 of the 
workflow process. 

The invention can be implemented or embedded in a 
variety of systems and devices to replace color interpolation 
25 technologies or subsystems. It can be included, for example, 
in implementation of the PostScript® language, or in a 
profile tool to populate a color lookup table in an ICC 
(International Color Consortium) device profile. An ICC 
profile is a description of the color characteristics for a 
30 device, such as a monitor, scanner or printer, and can be used 
to translate color data created on one device into the native 
color space of another device. To populate such a color 
lookup table, one needs color values at the grid points of the 
table. Generally, when taking measurements from color 
35 patches, one finds colors that are not on the grid of Lab 
colors, and therefore one has had to interpolate. With the 
present invention, one can calculate the values for a multi- 
dimensional color lookup table for an ICC profile and 
achieve the benefits pointed out in this specification. 
40 The invention can be implemented in digital electronic 
circuitry, or in computer hardware, firmware, software, or in 
combinations of them. Apparatus of the invention can be 
implemented in a computer program product tangibly 
embodied in a machine-readable storage device for execu- 
45 tion by a programmable processor; and method steps of the 
invention can be performed by a programmable processor 
executing a program of instructions to perform functions of 
the invention by operating on input data and generating 
output. The invention can be implemented advantageously 



appropriate colorants. In a PostScript® workflow, for 

example, the triple is taken as an index into the RenderTable 50 in one or more computer programs that are executable on a 



in a selected color rendering dictionary (CRD). The Ren- 
derTable is a three-channel color lookup table (CLUT). A 
linear interpolator calculates colorant values as a the result 
of a linear interpolation between two grid points of the 
RenderTable. In the case of an ICC (International Color 
Consortium) workflow, the triple is also used as an index 
into a CLUT (which is sometimes referred to as a calibration 
table) in an ICC device profile, and a linear interpolator 
calculates colorant values as a the result of a linear inter- 
polation between two grid points of the CLUT. 

As illustrated in FIG. 4, in a workflow in accordance with 
the present invention, a color lookup table 430 is used to 
produce an initial colorant vector 435 from an initial input 
color 420. The values at table grid points can be used, as can 



programmable system including at least one programmable 
processor coupled to receive data and instructions from, and 
to transmit data and instructions to, a data storage system, at 
least one input device, and at least one output device. Each 
55 computer program can be implemented in a high-level 
procedural or object-oriented programming language, or in 
assembly or machine language if desired; and in any case, 
the language can be a compiled or interpreted language. 
Suitable processors include, by way of example, both gen- 
60 eral and special purpose microprocessors. Generally, a pro- 
cessor will receive instructions and data from a read-only 
memory and/or a random access memory. Generally, a 
computer will include one or more mass storage devices for 



storing data files; such devices include magnetic disks, such 
values linearly interpolated from the grid point values, as are 65 as internal hard disks and removable disks; magneto-optical 
conventionally produced in the first part of a PostScript or disks; and optical disks. Storage devices suitable for tangi- 
ICC publishing workflow. Backward compatibility is pro- bly embodying computer program instructions and data 



06/29/2004, EAST version: 1.4.1 



US 6,646,763 Bl 

11 12 

include all forms of nonvolatile memory, including by way obtaining an initial set of colorant concentrations for the 

of example semiconductor memory devices, such as set of colorants, each concentration corresponding to a 

EPROM, EEPROM, and flash memory devices; magnetic respective colorant in a set of colorants; and 

disks such as internal hard disks and removable disks; generating iteratively a final set of colorant concentrations 

magneto-optical disks; and CD-ROM disks. Any of the 5 b y calculating a rendered color value, the rendered 

foregoing can be supplemented by or incorporated in, ^ yalue bd a ectral represe nt a tion of a 

ASICs (application-specific integrated circuits). reQdered coIof duced b a mbam& Qf colorants 

To provide for interaction with a user, the invention can . „ , , tl _ \ - . \ t . . . . 

. • i . . . . , . , denned by the set of colorant concentrations, calculat- 

be implemented on a computer system having a display , ' t t , , , ' , , 

device such as a monitor or LCD screen for displaying, in in S a difference between the rendered color and the 

information to the user and a keyboard and a pointing device 10 tar 8 et ^ a * d modifying the set of colorant concen- 

such as a mouse or a trackball by which the user can provide trations according to the calculated difference, 

input to the computer system. The computer system can be 2 - ^ method of claim 1, further comprising: 

programmed to provide a graphical user interface through receiving a threshold to be achieved, wherein the iteration 

which computer programs interact with users. terminates when a calculated difference between the 

By way of example, a printing device implementing an 35 rendered color and the target color is less than the 

interpreter for a page description language, such as the threshold or an iteration limit is reached. 

PostScript® language, includes a microprocessor for execut- 3. The method of claim 1, wherein obtaining an initial set 

ing program instructions (including font instructions) stored 0 f colorant concentrations comprises: 

on a printer random access memory (RAM) and a printer ^ the ta t color value t0 sdect an index int0 a 

read-ordymemory (ROM) and controlling a printer marking 20 colorant lookup table; and 

engine. The RAM is optionally supplemented by a mass . . . . . . . . , - . r . + 4 . f 

storage device such as a hard disk. Then essential elements obtamm g the , imtial ^ f ? ol ™ n ' concentrations from 

of a computer are a processor for executing instructions and A ^ e mde * ed contents of the lookup table, 

a memory. Acomputer can generally also receive programs 4 * ™ e method of claim wherem calculating a rendered 

and data from a storage medium such as an internal disk or 25 color corresponding to the set of colorant concentrations 

a removable disk. These elements will also be found in a comprises: 

conventional desktop or workstation computer as well as calculating a rendered color reflectivity function from the 

other computers suitable for executing computer programs set of colorant concentrations and a set of colorant 

implementing the methods described here, which can be reflectivity functions, each colorant reflectivity func- 

used in conjunction with any digital print engine or marking 30 tioD corresponding to a respective one of the colorants; 

engine, display monitor, or other raster output device anc * 

capable of producing color or gray scale pixels on paper, calculating a rendered color value from the rendered color 

film, display screen, or other output medium. reflectivity function. 

The invention has been described in terms of particular 5. The method of claim 4, wherein the rendered color 

embodiments. Other embodiments are within the scope of 35 value is calculated using an iiluminant selected in response 

the following claims. For example, the steps of the invention to a user input identifying desired viewing conditions, 

can be performed in a different order and still achieve 6. The method of claim 1, wherein modifying the set of 

desirable results. Functions other than 0 X ' that have the colorant concentration values comprises: 

property of being color additive can be used. Color- calculating a set of colorant concentration corrections 

difference formulas calculated in other spaces, such as 40 from the difference between the rendered color value 

CIELUV or C1E 94, with their corresponding partial and the target color value. 

derivatives, can be used. In addition, formulas developed for 7. The method of claim 1, wherein the target color value 

particular applications, such as the CMC(l:c) (Color Mea- is a triple defining a location in a color space, 

surement Committee) color difference formulas, with their 8. The method of claim 1, wherein the target color value 

corresponding partial derivatives, can be used. 45 defines a device-independent color. 

Device-independent color spaces other than the ones 9. The method of claim 1, wherein the difference between 

specifically mentioned can be used, such as calibrated RGB the rendered color and the target color is calculated in a 

and calibrated CMYK, each of which can be mapped to perceptually linear color space. 

XYZ triples and so can be accepted as input in place of XYZ 10. The method of claim 9, wherein the perceptually 
triples. Finally, while the invention has been designed for 50 linear color space is CIELAB or CIELUV. 
use in systems in which colorants are mixed physically, it 11. The method of claim 1, wherein the set of colorants is 
can also be applied where colorants are mixed visually, such an arbitrary set of two or more colorants, 
as by halftoning dots. Although the Kubelka-Munk theory 12. The method of claim 1, wherein the spectral repre- 
does not hold as well in the latter situation, other ink models sentation of each colorant derives from reflectivity measure- 
can be used for that situation; and in general the invention 55 ments of the colorant taken at intervals over the visible 
can be implemented on the basis of other ink models. spectrum without overprinting. 

In addition, the iterative process for refining a set of 13. The method of claim 1, wherein modifying the set of 

colorant vectors can be based on a technique for finding a colorant concentration values comprises: 

minimum of an error function other than one described modeling the difference between the target color and the 

above. For example, a Levenberg-Marquardt (LV) algorithm 60 rendered color as a scalar field over the set of colorants; 

can be used to solve the non-linear problem AE 2 =0. and 

What is claimed is: calculating a partial derivative of the surface with respect 

1. A method for color matching, comprising: to each colorant to calculate a correction for each 

receiving a target color value, the target color value being colorant. 

a nonspectral representation of a target color; 65 14. A computer program product, tangibly stored on a 

receiving a spectral definition of each colorant in a set of computer-readable medium, for matching a color, the prod- 
colorants; uct comprising instructions operable to cause a computer to: 



06/29/2004, EAST Version: 1.4.1 



US 6,646,763 Bl 



13 



14 



receive a target color value, the target color value being a 
nonspectral representation of the target color; 

obtain an initial set of colorant concentrations, each 
concentration in the set corresponding to a respective 
colorant in the set of colorants; and 

generate iteratively a final set of colorant concentrations 
by calculating a rendered color value, the rendered 
color value being a nonspectral representation of a 
rendered color produced by a mixture of the colorants 
defined by the set of colorant concentrations, calculat- 
ing a difference between the rendered color and the 
target color, and modifying the set of colorant concen- 
trations according to the calculated difference. 

15. The product of claim 14, wherein instructions to 
obtain an initial set of colorant concentrations comprise 
instructions to: 

use the target color value to select an index into a colorant 

concentration lookup table; and 
obtain the initial set of colorant concentrations from the 

indexed contents of the lookup table. 

16. The product of claim 14, wherein instructions to 
calculate a rendered color corresponding to the set of 
colorant concentrations comprise instructions to: 

select an illumination function; 



10 



20 



21. The method of claim 20, wherein the first set of 
colorant concentration values is derived by using the input 
color tuple as an index to obtain grid-point concentration 
values at two grid points of the color lookup table and 
calculating the first set of colorant concentration values as a 
linear interpolation of the grid-point concentration values. 

22. The method of claim 20, further comprising providing 
a color function table of color-additive function values, 
wherein the iterative non-linear process comprises: 

iteratively (a) calculating an interim color tuple from an 
interim set of colorant concentration values and the 
color function table, the initial interim set of colorant 
concentration values being the first set of colorant 
concentration values, and (b) deriving an interim set of 
colorant concentration values from a difference 
between the input color tuple and the interim color 
tuple. 

23. The method of claim 22, further comprising: 
calculating a partial derivative of an error function from 

the difference between the input color tuple and the 
interim color tuple; and 
using the partial derivative to derive a successor interim 
set of colorant concentration values from a current 
interim set of colorant concentration values. 

24. A system for mapping a nonspectral representation of 



calculate a rendered color reflectivity function from the 25 an input color to a set of concentration values for a set of 



set of colorant concentrations and a set of colorant 
reflectivity functions, each corresponding to a respec- 
tive one of the colorants; and 
calculate a rendered color value from the illumination 
function and the rendered color reflectivity function. 

17. The product of claim 16, further comprising instruc- 
tions to select the illumination function in response to a user 
input identifying a desired viewing condition. 

18. The product of claim 14, wherein instructions to 
modify the set of colorant concentration values comprise 
instructions to: 

calculate a set of colorant concentration corrections from 
the difference between the rendered color value and the 
target color value; and 

add the set of colorant concentration corrections to the set 
of colorant concentrations. 

19. A system, comprising: 

means for receiving a target color value, the target color 
value being a nonspectral representation of a target 
color; 

means for obtaining an initial set of colorant 
concentrations, each concentration in the set corre- 
sponding to a respective colorant in the set of colorants; 
and 

means for generating iteratively a final set of colorant 
concentrations by calculating a rendered color value, 
the rendered color value being a nonspectral represen- 
tation of a rendered color produced by a mixture of the 



35 



45 



50 



device-specific colorants, the representation being an input 
color tuple, comprising: 

means for using the input color tuple to derive a first set 

of colorant concentration values from a color lookup 

table; and 

means for refining the first set of colorant concentration 
values by an iterative non-linear process to generate a 
final set of colorant concentration values. 

25. The system of claim 24, wherein: 

the system further comprises means for providing a color 
function table of color-additive function value; and 

the iterative non-linear process comprises means for itera- 
tively performing the following actions: (a) calculating 
an interim color tuple from an interim set of colorant 
concentration values and the color function table, the 
initial interim set of colorant concentration values 
being the first set of colorant concentration values, and 
(b) deriving an interim set of colorant concentration 
values from a difference between the input color tuple 
and the interim color tuple. 

26. The system of claim 25, further comprising: 
means for calculating a partial derivative of an error 

function from the difference between the input color 
tuple and the interim color tuple; and 
means for using the partial derivative to derive a succes- 
sor interim set of colorant concentration values from a 
current interim set of colorant concentration values. 

27. A computer program product, tangibly stored on a 



colorants defined by the set of colorant concentrations, 55 computer-readable medium, for mapping a nonspectral rep- 



calculating a difference between the rendered color and 
the target color, and modifying the set of colorant 
concentrations according to the calculated difference. 
20. A method for mapping a nonspectral representation of 
an input color to a set of concentration values for a set of 
device-specific colorants, the representation being an input 
color tuple, comprising: 

using the input color tuple to derive a first set of colorant 

concentration values from a color lookup table; and 
refining the first set of colorant concentration values by an 
iterative non-linear process to generate a final set of 
colorant concentration values. 



60 



65 



resentation of an input color to a set of concentration values 
for a set of device-specific colorants, the representation 
being an input color tuple, the product comprising instruc- 
tions operable to cause a computer to: 
use the input color tuple to derive a first set of colorant 

concentration values from a color lookup table; and 
refine the first set of colorant concentration values by an 
iterative non-linear process to generate a final set of 
colorant concentration values. 
28, The product of claim 27, further comprising instruc- 
tions operable to cause a computer to derive the first set of 
colorant concentration values by using the input color tuple 
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as an index to obtain grid-point concentration values at two 
grid points of the color lookup table and calculate the first set 
of colorant concentration values as a linear interpolation; of 
the grid-point concentration values. 

29. The product of claim 27, for use with a color function 
table of color-additive function values, the product further 
comprising instructions operable to cause a computer to 
perform the iterative non-linear process iteratively (a) cal- 
culating an interim color tuple from an interim set of 
colorant concentration values and the color function table, 
the initial interim set of colorant concentration values being 
the first set of colorant concentration values, and (b) deriving 
an interim set of colorant concentration values from a 
difference between the input color tuple and the interim 
color tuple. 
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30. The product of claim 29, further comprising instruc- 
tions operable to cause a computer to: 

calculate a partial derivative of an error function from the 
difference between the input color tuple and the interim 
color tuple; and 

use the partial derivative to derive a successor interim set 
of colorant concentration values from a current interim 
set of colorant concentration values. 

31. The product of claim 30, wherein the input color tuple 
is a device-independent triple and the interim color tuple is 
a device -independent triple. 
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